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Abstract 

A completely integrable dynamical system in discrete time is stud- 
ied by means of algebraic geometry. The system is associated with 
, factorization of a linear operator acting in a direct sum of three linear 

' spaces into a product of three operators, each acting nontrivially only 

in a direct sum of two spaces, and the following reversing of the or- 
der of factors. There exists a reduction of the system interpreted as a 
classical field theory in 2 + 1-dimensional space-time, the integrals of 
motion coinciding, in essence, with the statistical sum of an inhomoge- 
^\ , neous 6-vertex free-fermion model on the 2-dimensional kagome lattice 

(here the statistical sum is a function of two parameters). Thus, a 
connection with the "local" , or "generalized" , quantum Yang — Baxter 
O 1 equation is revealed. 

D ! 

^ c| This paper is devoted to an exactly solvable model in discrete time asso- 

ciated with factorization of a certain matrix in a product of three matrices 
and the following change of the order of those matrices into the inverse one, 
and also to a reduced case of this model, when it becomes a sort of field the- 
ory in 2 + 1-dimensional wholly discrete space-time. The generating function 
of the motion integrals in this reduced case is nothing else but a statistical 
sum of an inhomogeneous 6-vertex model (depending on 2 parameters) on 
the kagome lattice. One can say thus that the reduced dynamical system 
describes an evolution in discrete time of the 6-vertex model Boltzmann 
weights which conserves the statistical sum. The evolution is local, which 
means that the value of each weight in the moment r + 1 is affected only 
by (itself and) several neighboring weights in the moment r. The weights 
satisfy the "free fermion" condition known from the quantum Yang — Baxter 
equation (QYBE) theory. 

The evolution of weights can be seen as a certain realization of ideas of 
the paper [B| concerning the "localization" of such equations as QYBE. As 
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is known, the usual QYBE 



R12L13M23 — M23L13R12 

contains the same operators in the L.H.S. as in the R.H.S. (the indices 
pointing to different linear spaces in whose tensor product an operator acts). 
The local QYBE 

R 12 Li 3 M 2 3 = M 23 L l3 R l2 (1) 

has different operators in each side. Note, however, that the fundamental 
equation in the present work is an equation differing from (j^) in that the 
operators in it act not in a tensor product, but in a direct sum of three 
spaces! This will be seen in Section [l], where we will have to represent 
as a product of three matrices some matrix obtained at the previous step 
of evolution as a similar product, but in the inverse order. How a tensor 
product arises from a direct sum, can be explained by using the Clifford 
algebras. However, in this paper it is not done, and the 6-vertex model is 
obtained through some consideration of topological character. 

Note also that the connection between the statistical physics models of 
free-fermion type and completely integrable models is known for a long time 
now — see and the works cited therein. The first of the author's works 
on this subject was 110], where a statistical sum conserving evolution of the 



weights was proposed for a flat dimer model. In both papers, [10[ and the 
present one, the tool of investigation is algebro-geometrical objects: alge- 
braic curves and divisors. These objects are well known in mathematical 
physics since the works on finite-gap integration of soliton equations ap- 
peared (see e.g. the book [||). In dealing with them, I try to hold to a 
reasonable degree of strictness, in order not to overload the text with un- 
necessary details. 



1 Definition of dynamical system 

Let 

A = 




be a block matrix acting in the linear space of m+n+r-dimensional complex 
column vectors, so that, e.g., A, F and K are square matrices of sizes m x 
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m, nxn, r xr respectively. Consider the problem of factorization of matrix 
A into a product of the form 



M 







Ci 













1 




^3 


-( 




A = ^1^2-43, (3) 

where 

-4j - I Ci Dj 1 . A-2 

I Co Do I 

(4) 



A 3 B 3 
C 3 D 3 

The factorization (||) may be seen as a generalization of the factorization 
of an orthogonal rotation in the 3-dimensional space into rotations through 
the "Euler angles" . However, no orthogonality conditions will be considered 
in this paper. 

One can obtain from the factorization (^|) other factorizations of the 
same kind by using the following transformation of the triple Ai, A2, A 3 : 



Ai^Ai\ M 2 I , ^3 -»■ I Ah' 1 I A 3 
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(5) 

i 00\ / 10 °\ 

A-2 - I 


Mi, M2 and M3 being arbitrary non-degenerate matrices of proper sizes. 

Lemma 1.1 For a generic matrix A factorization (^) is (if exists) unique 
to within the transformations 

Proof. First, let us specify that matrix A will be called generic with 
respect to this lemma if it a) is non-degenerate and b) allows a factoriza- 
tion @ with matrices A2 and V2 non-degenerate. Let 

A = A'^A'z (6) 

be another factorization of the same kind. From @ and (^) one finds 

a' 2 = a';a 2 ai (7) 
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where 



A'{ = (AW'Au A = A 3 (A' 3 ] 



(8) 



Let us denote the blocks in the dashed matrices by the same letters A\, . . . , 
-D3 as in equalities (Q) with proper number of dashes added to them. The 
relation (0) is rewritten as 



A' 2 B' 2 
10 

C 2 D' 2 



( A'{A 2 A'{B 2 C'{ + B'lA'i A'{B 2 D' 3 ' + B'{B'i 
C'{A 2 C'{B 2 C'i + D'(A'i ClB 2 D'i + D'{B% 
\ C 2 D 2 C'i D 2 D'i 



(9) 



From here, one obtains at once the equalities 

C 2 = C 2 , C[ = 0, C 3 



0. 



Taking this into account, one finds from the block in 2nd row and 2nd 
column that 

D'{A'i = 1. 

Thus, D'l and A'l are non-degenerate. Now the blocks just above the main 
diagonal yield 

0, B'i = 0. 



B>> 



So, the matrices A" and A'% (g) are block-diagonal. In is easy to see that 
this means exactly that A[ , A 2 , A' 3 are obtained from A\ , A 2 , A3 by the 
transformation (|5|). The lemma is proved. 

Now let us construct, starting from the block martix A, new matrix B 
by following means: factorize A into the product (||) and set 



B = A 3 A 2 Ai. 



(10) 



From the above considerations it is seen that the matrix B is determined to 
within the transformations 



Let us call such transformations, as applied to the block matrices here, the 
gauge transformations. The following simple but important observation is 
valid: if matrix A itself undergoes a gauge transformation, this in no way 
affects the set of matrices B obtained from formula (|7^). 
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M 2 - 1 
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M 3 - 
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It will be shown in Section || that factorization (|3|) does exist for a generic 
matrix A. This factorization will be constructed by means of algebraic ge- 
ometry. Taking this into account, we are ready now to define the dynamic 
system that we are going to examine. Let A4 be the set of block matri- 



ces (|J) taken to within gauge transformations (11), or, using stricter lan- 
guage, the set of equivalence classes of such matrices with respect to trans- 
formations (|TT|) . The set Ai will be our "phase space". Then, the birational 
mapping / is defined on the set M that brings into correspondence with 
a matrix .4, factorized into the product (0), the matrix B factorized into 
the product (p7o|) . Let us now bring into consideration the "discrete time" r 
taking integer values and say that to the trasition from time r to time r + k 
corresponds the mapping 

k times 

2 Invariant algebraic curve of matrix A and some 
divisors on it 

The dynamical system of Section |l] turns out to be completely integrable. To 
be exact, an invariant curve T can be constructed out of matrix A, together 
with a divisor T> on it. In terms of these algebro-geometrical objects, the 
evolution is as follows: T does not change, while T> — more precisely, its linear 
equivalence class — depends linearly on the discrete time r. 

Let us start from the definition of the curve T. The word "invariant" in 
this definition will be justified in Section ||. 

Definition 2.1 The invariant curve T of the operator A of the form ^) is 
an algebraic curve in CP 1 x CP 1 x CP 1 (i.e. in the space of three complex 
variables u,v,w, each allowed also to take value oo) given by equations 



det(A 



( ul m 

vl n |) = 0, (12) 




V 

v = uw. (13) 



Here a subscript of each 1 means the size of corresponding unity martix, 
while denotes rectangular zero matrices of different sizes. Strictly speak- 
ing, equations (|l^, |l|) define the "finite part" of the curve T, the whole 
curve T being its closure in Zariski topology. 
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( 
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uX \ 


A 








vY 




\ 




\ 


wZ J 



X 

The equality (|l^) means that a column vector X = | Y | exists, with 

Z 

X, Y, Z column vectors of dimensions m, n and r correspondingly, such that 



(14) 



Such vectors X form a one-dimensional holomorfic bundle over T. 

The next lemma shows the structure of the zero and pole divisors of 
functions u,v,w. For these divisors, the notations (it)o, (u)^ etc. are used. 

Lemma 2.1 There exist such effective divisors (i.e. finite sets of points) 
T>i , . . . , T>q in the curve T that 

(u) 00 = V 1 + V 2 , (v) 00 = V 1 + V 3 , {w) 00 =V- i + V A , 

{u) = V i + V 6 , (v) = V 5 + V e , (w)o = V 5 + V 2 . [ ' 

T>3 and P5 are of degree m, T>2 and T>^ are of degree n, T>\ and T>§ are of 
degree r. Generally, all points included in divisors T>i, . . . ,T>q are different 
from each oter. 

Proof. Consider, e.g., the case u = 0, w ^ 0, w ^ 00. Then, according 
to (|l3|), v = 0. The equality ( |14| ) turns into the following system: 

AX + BY + CZ = 0, 
DX + FY + GZ = 0, 
HX + JY + KZ = wZ. 

One can express X and Y through Z (e.g., Y = —(F — DA^ 1 B)~ 1 (G — 
DA~ 1 C)Z) and then substitute these expressions into the third one. One will 
come to an equation of the form 

KZ = wZ (16) 

which has r characretistic roots w\, . . . , w r , different from each other in 
general case. This is how r points (0, 0, w\), . . . , (0, 0, w r ) of divisor Vq are 
obtained. The other divisors in ([l^) arise in a similar way. The lemma is 
proved. 
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( x 

The vector Y | in ([lj) is determined up to a scalar factor which 

v z 

may depend on the point in the curve T. So, this vector can be normalized 

( X 

by setting its first coordinate identically equal to unity (cf. Q). Y 

V z 

becomes then a meromorphic vector on V with a certain pole divisor T>. 
However, X, Y and Z taken separetely satisfy stronger restrictions, as the 
following lemma shows. In the lemma, (/) denotes the divisor of a function 
/ (zeros enter with the + sign, poles with the — sign, as usual). 

Lemma 2.2 The column vector X consists of functions f such that 

(/)+ 2? -(«)«, >0; (17) 

the column vector Y consists of functions f such that 

(f)+V-(v) oo >0; (18) 

the column vector Z consisits of functions f such that 

(f)+V-(w) ao >0. (19) 

Proof. One can see immediately from the formula ( |i~4| ) that the vector uX 

( X ' 

entering into R.H.S. cannot grow faster than the vector 



Y I in L.H.S. in 

Z 

such points where u = oo. This is exactly what the inequality (|17|) states. 
The inequalities (18) and (|l9|) are proved similarly. 



V 



Now the time has come to make it sure that the curve T, for a generic 
matrix A, is a smooth irreducible curve. One may wish also to calculate 
its genus in some simple way. To do that, we are now going to examine 
a relatively simple particular case of the matrix A, although at the time 
"generic" enough to make sure that such its features as genus and the degree 
of divisors are the same for matrices in some Zariski neighborhood. 

Thus, let all the matrix elements of A equal zero except the ones lying, 
first, in the main diagonal and, second, in the "broken" diagonal parallel to 
the main one (for these latter matrix elements, the difference between the 
numbers of a column and a row must be some constant modulo m + n + r). 
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The elements in the main diagonal will be denoted as a%, . . . a m , fx, . . . f n , 
kx,...k r ; and let the elements in the broken diagonal be all equal to the 
same complex number s: 

( «i 



A 




\ 



(20) 



k r J 



It does not matter through which blocks exactly the "broken" diagonal 
passes. 

For the finite u, w, the curve T now examined is given by equation (re- 
sulting from the substitution of into (|l2j)) 



in 



n 



F(uw) = Yl 
As is known, in singular points 



a=l 0=1 



uw) ■ J] {ky - w) ± s m+n+r = 0. (21) 

7=1 



dF/du = 0, 
dF/dw = 0. 



(22) 



The system (22) has a finite number of solutions, and, changing s in (21), 
one can make these solutions not to lie in the curve T, which thus will be 
free of singularities for finite u, v . It is an easy exercise to show that there 
are no singularities when u or w is infinite as well. 

Returning now to general martices A and curves T, let us note that 
it cannot be that the system p2] ) or its substitute in the neighborhood of 
infinite u or w possesses solutions in the curve T in general case and does not 
possess them in a particular case. Thus, the smoothness of T for a generic A 
is clear. As for irreducibility, to prove it let us examine the natural projection 
of r onto the Rimann sphere CP 1 of the variable u. This projection is an 
n + r-sheet cover, and if V consisted of two or more components, the sheets 



8 



of the cover would split into groups belonging to each component. To prove 
that it is not so in the general case, it is enough to present an example 
where it is not so. To do this, take A of the form ( p0| ) and, moreover, put 
fx = . . . = f n = k\ = . . . = k r = 0. Equation (^T|) then becomes 



rn 



w n+r u n 



n k - u) ± s m+n+r = o. 



a=l 



Let ax ^ and not coinside with other a a . Then in a neighborhood of the 
point (u,w) = (ai,oo) w behaves like 



w 



1 ~ {ax - u) l ^ n+r \ 



From here one sees that all the mentioned n + r sheets belong to a single 
component, i.e. the irreducibility of T is proved. 

Now let us denote the number of branch points of the cover r — > C 3 u 
as b. Then the genus of the curve T, according to the Riemann — Hurwitz 
formula, is 

g = \-n-r+ h -. (23) 
Our next aim is to express b and g through m, n and r. 



Lemma 2.3 The degree of the vector 



( X 

Y | pole divisor T> is m + 6/2. 

V 



z 

Proof. Write out the equation ( |T4| ) "explicitly": 
AX + BY + CZ = uX, 

DX + FY + GZ = vY, } (24) 
HX + JY + KZ = wZ. 

Expressing X through Y and Z by means of the first of these equations and 
substituting into the rest, one finds: 

( u- 1 (D{u- A)- 1 B + F) u- 1 (D{u- A)- l C + G)\( Y\_ (y\ 
\ H{u-A)- l B + J H{u-A)' l C + K )\ Z ) ~ W \ Z J ' 

(25) 

To a generic u correspond n + r different w = wx, ■ ■ ■ w n + r , and the corre- 

( Y \ 

sponding n + r vectors I J are linearly independent as eigenvectors of 
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the matrix in L.H.S. of (p5|). An easy check shows that in the "suspicious", 
from the standpoint of equation (|25|), points u = 0, oo, and such points 

' Y \ 



where det(u — A) = 0, there exist n + r linearly independent vectors 



Z 



as well. 

Consider a determinant 



d 



Y(wi) ... Y(w n+r ) 
Z(wi) . . . Z(w n+r ) 



Given u, it changes its sign under odd permutations of w's. This means 
that d 2 is a function of u only. From the above one sees that the number of 
zeros of function d 2 {u) equals b, because it is in the branch points and only 
in them that d 2 (u) vanishes. Thus, the number of function d 2 (u) poles also 

( Y \ 

equals 6, and the pole divisor degree of the meromorphic vector I J is 

6/2. All this consideration is perfectly standard, see [Q, g|. 

Finally, from the 2nd and 3rd equations of the system ( p4j ) one sees that 
the vector X has its poles, when Y and Z are finite, in those points of T 



where v = w = oo and only in them. Formulae (15) show that these are the 
points of divisor P3, and there are m of them. This adds m more poles to 
6/2 already present, and with this the proof of the lemma is complete. 
Let us turn again to matrices A of the form (|20|). For such matrix, it is 

( x \ 

easy to find the vector Y = X in a given point (u, v, w) G T. Let us 

V z J 

assume that the "broken" diagonal is the one adjacent to the main diagonal, 
so that there is only one letter s in the lower left corner. Then the following 
holds for the vector X coordinates: 

(01 - u)X x + sX 2 = 0, 
(o 2 - u)X 2 + sX 3 = 0, 



(dm ~ 


- u)X m + sYi = 0, 


(fx- 


v)Yi + sY 2 = 0, 


(fn- 


v)Y n + sZx = 0, 


(h- 


w)Zi + sZ 2 = 0, 



(k n - w)Z n + sXi = 0. 
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From here the ratios between vector X coordinates are readily seen. 
Assuming the normalization condition X\ = 1, one finds out that X has 
the poles a) of the order ra in n points (u,v,w) = (oo,fp,0) and b) of 
the order ra + n in r points (u,v,w) = (00,00, & 7 ). In all, X possesses 
thus mn + mr + nr poles, taking their multiplicities into account. Recalling 
Lemma |2.3| and formula (|2^ ) , one can now find the genus g of the curve as 
well. As a matrix A of the form (|20| ) is "generic enough" , the results on the 
degree of divisor T> of the vector X and genus g of the curve apply also to 
curves corresponding to generic matrices A. Let us formulate them as the 
following lemma. 



Lemma 2.4 For a generic matrix A, the genus of the curve T is 
g = ran + mr + nr — ra — n — r + 1, 



while the degree of divisor V of the meromorphic vector X 



( X 
Y 
Z 



(26) 



is 



ran + rar + nr = g + ra + n + r — 1. 



(27) 



Thus, in this section we have constructed, for a given matrix A, an 
algebraic curve T and a bundle of vectors X over it, and calculated the 
genus g of the curve and the degree of the bundle (i.e. the divisor T> degree). 
As a helpful tool, a matrix A of special simple form (2C) was used which, 
from many viewpoints, was "generic enough". In Section || we will study 
how these objects behave under evolution introduced in Section |l|. 



3 Evolution in terms of divisors 

In this section it is shown, at first, that there exists a one-to-one correspon- 
dence (more precisely, a birational isomorphism) between the set of block 
matrices A (||) taken up to gauge transformations @ , and the set of pairs 
(an algebraic curve, a linear equivalence class of divisors on it) of a certain 
kind. This correspondence has, in essence, been constructed in Section 
and here are some missing details. Then, it is explained which divisors and 
why correspond to the factors A\,A2 and ^3 in (H) taken separately. Fi- 
nally, it is demonstrated that to the matrix B (|l0|) obtained from A by 
reversing the order of its factors, the same curve V corresponds, but the 
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divisor undergoes some constant shift. Thus, the motion linearizes in the 
Jacobian of curve T. Let us proceed to a detailed consideration. 

Equations ( |l2| , |T3| ) define, for a block matrix A, an algebraic curve T. 
Those equations can obviously be written as 



m n r 



J2J2J2 CL ijkU i V j W k = 0, 
1=0 j=0 k=0 

(28) 

V = uw. 

Besides, a linear bundle over T has been constructed in Section — the bun- 

( x \ 

die of vectors Y (0). That means that the divisor T> of the bundle 

is determined, up to linear equivalence, whose degree is g+m+n+r— 1, 
g being the curve's genus (27). Gauge transformations (||) do not change a 



pair (r, class of divisor T> ). 

Now let us show how to construct the matrix A starting from coefficients 
ajjfc of the curve (^) (arbitrary complex numbers in "general position" ) and 
a divisor T> of degree g + m + n + r — 1. Note that genus g of the curve 



r defined by formulae (25) without any (a priori) connection with block 
matrices is given by the same formula (|26|). This can be seen, e.g., by 
starting again from the "simple" curve (|2l]) of Section |2| whose genus is 
known. Define now the meromorphic column vectors X, Y and Z, guided 
by Lemma for components of vector X, take m linearly independent 
meromorphic functions on V satisfying relation (|l7|), and for Y and Z take, 
similarly, n functions satisfying (|l^) and r functions satisfying (|i~9|). 



Note also that Lemma 2.1 about divisors (n)oo, (u)oo, {w)oo entering in 
formulae (17-1^) remains valid for curves defined by an "abstract" system 
(^8|), which is immediately seen on substituting zero or infinity for u,v, or 
w in d|). 

It is clear now that relation (|14|) determines unambiguously the matrix 
A (cf . a similar construction in paper Q ) . Another choice of linearly inde- 
pendent functions for components of X, Y and Z leads, of course, to a gauge 
transformation <M\. 



Examine now each multiplier in factorization (f|) separately. Lemma 1.1 
shows that factorization (||), if exists, is unique to within the transforma- 
tions (||) . Let us demonstrate how to construct this factorization by algebro- 
geometrical means. 
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Figure 1: Factorization of matrix A and the divisors 
Consider the following figure (Fig. |]). 

The meaning of the numbers standing near the edges in this figure is as 
follows: if those numbers are jk, then the meromorphic vector corresponding 
to the edge consists of such functions / whose zero and pole divisor (/) 
satisfies inequality 

(/) + V - V, -V k >0. 



Those inequalities must be in agreement with Lemmas and 2.2. In par- 
ticular, the matrix ^3 will be defined by equality (notations of formulae (|j) 
are used), 

(££)(!;)-(£)■ 

where the meromorphic vector Y consists of functions / such that 

(/) + V - £>i - £> 3 > 
(formulae (|l~8|) and (|l~5|)); Z of functions such that 

(f)+V-V 3 -Vt >0 
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(formulae ( |l9D and fll5|)); Y' and Z' consist by definition of such linearly 
independent functions that 

(f)+V-V 3 -V 6 >0 

for Y' , and 

(f)+V-V 2 -V 3 >0 
for Z'. It is easy to see that ( p9| ) is a correct definition for matrix .A3, 

because the components of each of the vectors ^ ^ J and ( ^ ) form a 

basis in the space of meromorphic functions / such that 




(/) + x> - x> 3 > 0. 



Next, let 





(30) 

Q x Vl 11 V I 1 ,.V I ' ( 31 ) 

where X' consists of functions / such that 

(f) + V-V 2 -V e > 0. 

It is shown in much the same way as above that equalities ( |30| ) and ( |3l| ) 
do correctly define the matrices A2 and A\. What remains is to check the 
validity of equality (||) for A\, A2 and ^,3 given by these definitions. To do 
this, observe that (29-|31~D together yield 







( uX 




- 


vY 






{ wZ 



A!A 2 A 3 I Y J = I vY J . (32) 

The equality (^) follows from comparing (p^) with (p^). 

Note that the arbitrariness in choosing X',Y' and Z' corresponds, of 
course, to transformations (|5|). 

Now let us pass to matrix £>, a product of the same three factors in the 
inverse order. The formulae ([H]) and ( p9| - |3"T|) t ogether yield (if one multiplies 
both sides of ( |30| ) by u, and both sides of (|29|) by v ): 



x' ^ 






Y' 


I" 




uZ' 1 







iB y = «y . (33) 



14 



Compare the divisors of meromorphic vectors in L.H.S.'s of (|33j) and (14). 
An easy calculation shows that 



V x . -V x = V Y , -Vy = V (uZI) -Vz = V x - V 6 . 

One sees hence that the same curve T corresponds to the operator B as to 
the operator A, while the divisor T> changes to T> + T>\ — Vq. 

Thus, in this section the name "invariant" has been justified for the curve 
T: it has been shown not to change under the evolution of Section |]. At 
the same time, it was demonstrated how to construct the factorization (H). 
Finally, it was shown that the evolution is described in algebro-geometrical 
terms as a linear, with respect to discrete time, change of (the linear equiv- 
alence class of) divisor V: it changes by V\ — Vq per each unit of time. 

4 Connection with the 6-vertex model on the 
kagome lattice 

Consider now a reduction of the system defined in Section [l| which leads 
to a dynamical system in 2 + 1-dimensional fully discrete space-time. Let the 
linear space in which matrix A @ acts have a basis enumerated by edges of 
a triangular lattice on the torus (Fig. |2|) . Let the lattice contain m horisontal 
edges, n oblique edges and r vertical edges, and let this correspond to the 
block structure (||) of matrix A in the sense that, say, the block A acts in 
the space generated by basis vectors corresponding to horisontal edges and 
so on. Impose the following "locality" condition on matrix A: let the vector 
corresponding to any given edge of the lattice be transformed under the 
action of A into a linear combination of just three vectors, corresponding to 
the edges coming upwards, to the right and northeastwards from the vertex 
that is the upper, right, or northeastern end of the considered "incoming" 
edge (Fig. |3|). Thus, only those elements of matrix A are not zeros that 
correspond to "local" transitions of Fig. |3[ 

The factorization of a "local" matrix A into the product (||) corresponds 
to each vertex represented by a small circle in Fig. [2] being converted into 
a triangle of the type shown in Fig.|l], so that the lattice transforms into a 
kagome lattice (Fig. ||). The triangles arising from the vertices-circles are 
shaded in Fig. |I[ 

A step of the evolution consists in those triangles being "turned inside 
out" , or else one may say that each oblique line replaces the one above it. We 



15 



Figure 2: The triangular lattice 




Figure 3: The vectors corresponding to "incoming" edges are transformed 
by A into linear combinations of those corresponding to "outgoing" edges 




Figure 4: The kagome lattice 
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w^ 1 



will see soon that this procedure is connected with solutions of the "local" 
in the sense of || quantum Yang — Baxter equation. 

Now return to the triangle lattice of Figure |2[ Let us express the "integral 
of motion" 



(34) 



in terms of paths going along the edges of this lattice (I(u, w) is indeed an 
integral of motion with any u, w, because the equality I(u, w) = determines 
the invariant curve, and a possible multiplicative constant is fixed by the 
fact that the constant term in ( |34"| ) equals unity). 

As it known, the determinant of a martix is an alternating sum of its 
elements' products, each summand corresponding to some permutation of 
the matrix columns, while each permutation factorizes into a product of 



cyclic ones. As applied to our matrix A, it means that the determinant (34) 
is a sum each term of which corresponds to a set of closed trajectories going 
along the arrows according to Fig. |3] (recall that the lattice is situated on 
the torus!). The trajectories of each given set can have intersections and 
self-intersections, but none of the edges may be passed through twice or 
more by one or several trajectories. 

To be exact, to each trajectory corresponds a product of entries of the 
matrix 





( 


u- 1 





A 











\ 









w 





-i 



corresponding to transitions through a vertex to a neighboring edge accord- 
ing to Fig. ||, multiplied (the product as a whole) by (—1). To each set of 
trajectories (including, of course, the empty set) corresponds the product 
of the mentioned values corresponding to its trajectories. The reader can 
verify that all the minus signs, including that in formula (|34"|), are taken into 
account correctly. 

It is easy also to describe the determinant I(u, w) in terms of the kagome 
lattice obtained on factorizing the matrix A into the product (0). This 
description almost repeats two preceding paragraphs. Let us formulate it as 
the following lemma. 

Lemma 4.1 I(u,w) is a sum over sets of trajectories on the kagome lat- 
tice; the direction of motion is upwards, to the right, or northeastwards; 
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none of the edges is passed through twice by trajectories of a given set; to 

the vertices of types 7^, I , X , if a trajectory passes through them, 
correspond the factors equalling matrix elements of matrices Ai, A2, A3 
respectively; besides, to each move to the right through a lattice period corre- 
sponds a factor u~ 1 , and to each move upwards — a factor w^ 1 (and both of 
them to a diagonal move); finally, to each set corresponds one more factor, 
f—\\ (number of trajectories) 

Now let us link I(u, w) with the statistical sum of (inhomogeneous) 6- 
vertex model on the kagome lattice. Let each edge of the kagome lattice 
be able to take one of two states, which will be depicted below as either 
presence or absence of an arrow on the edge (the arrow will always be di- 
rected upwards, to the right, or northeastwards). A "Boltzmann weight" 
will correspond to each vertex as follows: if there are no arrows on the edges 
meeting at the vertex, the weight will be 1; if there is exactly one arrow 
coming into the vertex and exactly one going out of it, the weight will be 
equal to the corresponding matrix element of A\, A2 or .4.3 (e.g. to the 



vertex 



corresponds the weight equal to the matrix element of A2 



that is responsible for a transition between the vector corresponding to the 
left edge, and the vector correcponding to the upper edge); if there are 2 
incoming and 2 outgoing arrows, the weight is the difference between the 
products of weights corresponding to the intersecting and non-intersecting 
paths through the vertex: 



Weight 



Weight 



Weight 



-Weight 



Weight 



(35) 



in the rest of cases the weight is zero. 

A weight will also correspond to each edge of the kagome lattice: weight 
1 to an edge without an arrow, and weights m" 1 / 2 , w~ 1 / 2 or u~ 1 / 2 w~ 1 / 2 to 
a horizontal, vertical or oblique edge having an arrow. If needed, the edge 
weights can be included in the vertex weights, but we will not do that here. 

The statistical sum S(u,w) of our 6-vertex model is, of course, a sum of 
products of vertex and edge weights over all arrow cinfigurations. The next 
lemma is the key statement. 
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Lemma 4.2 The statistical sum S(u, w) is a sum over the same sets of tra- 
jectories as the determinant I(u,w), and to each set corresponds the same 
summand up to, maybe, a minus sign. To be exact, the number of tra- 
jectories in the exponent of ( — 1) in Lemma \4-l\ changes to the number of 
intersections (self-intersections included) of a given set of trajectories. 

Proof is evident from the statistical sum definition. 

Each closed path on the torus is homologically equivalent to a linear 
combination of two basis cycles a and b. The same is true for a set of 
paths (trajectories), regarded as a formal sum of them. Different sets may 
be homologically equivalent to a given cycle la + mh but, as the following 
lemma shows, they have something in common. 

Lemma 4.3 For any set of trajectories on the torus homologically equiva- 
lent to a cycle la + mb (a, b being basis cycles, l,m — integers), 

(number of intersections) — (number of trajectories) = Im — I — m{mod2). 

(36) 

Proof may consist in the following simple consideration: 1) if the set 
consists of I trajectories going along a, and m ones going along b, (^) is 
obviously true, 2) under deformations of trajectories, the number of intersec- 
tions changes only by even numbers, 3) with elimination of an intersection 




or inverse operation, the L.H.S. of (|36|) may change 



also only by an even number. Starting from an arbitrary set and applying 
the transformations 2) and 3), one can arrive at a set of type 1), so the 
lemma is proved. 

Let the following products of edge weights correspond to the basis cycles: 
x = u ai w^ 1 for a and y = u a2 w@ 2 for b. Denote 

s(x, y) = S(u, w), f(x, y) = I(u, w). (37) 



Theorem 4.1 The statistical sum of the inhomogeneous 6-vertex model on 
the kagome lattice defined in this section is invariant with respect to the evo- 
lution of the reduced 2 + 1-dimensional model (for all u,w) and is connected 
with the determinant I(u, w) §34\), whose vanishing defines the invariant 
curve of the model, by relations ( in the notations of ( p^ j ) 

s(x, y) = 1/2 {-f{x, y) + f(-x, y) + f(x, -y) + f(-x, -y)) , (38) 

fix, y) = 1/2 (s(x, y) + s(-x, y) + s(x, -y) + s(-x, -y)) . (39) 
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Proof. It follows from Lemmas 4.1, 4.2, |4.3| that in the expansions of 
s(x,y) and f(x,y) in powers of x and y the coefficients near x l y m coincide 
if Im — I — m is even, and differ in their signs in the opposite case. This is 
exactly what the formulae (35,39) are about. 



5 Discussion 



The reader can verify that the 6-vertex model weight matrices correspond- 
ing to vertices of a triangle of the kagome lattice, and those matrices cor- 
responding to the same triangle "turned inside out" , indeed satisfy a local 
Yang — Baxter equation of the type of (||) . Such local, or generalized^ equa- 
tions are gaining popularity, and now not only solutions to the generalized 
Yang — Baxter equation, but even those for generalized tetrahedron equation 
have been found j|, ||, § . The author hopes to show soon that the construc- 
tions in this paper do have relationship with the tetrahedron equation and 
its generalizations. 

It is interesting also to link our model on the kagome lattice with the 
model generated by the star — triangle transformation ||] and described in 
the end of the paper |pX| ]. It turns out that, for that purpose, one has at first 
to examine the case of orthogonal matrix A (§J3|) (its factors A\ , A2 , -4.3 being 
orthogonal, too). The orthogonality corresponds to interesting symmetry 
conditions imposed on the invariant curve V and divisor T>: the canonical 
class of r is needed to formulate the condition on T> in the most elegant 
way! The details about the star — triangle model, orthogonality etc. will be 
presented in one of the author's coming publications. 

Thus, it was a certain generalization of the star — triangle transforma- 
tion that has been presented in this paper, and one more aspect of close 
connection between the classical and quantum integrable models has been 
demonstrated. 
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